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Abstract 

Extended theories are widely used in the literature to describe the relativistic fluid. The mo- 
tivation for this is mostly due to the causality issues allegedly present in the first order theories. 
However, the decay of fluctuations in the system is also at stake when first order theories that 
couple heat with acceleration are used. In this paper it is shown how the generic instabilities in 
relativistic fluids are not present when a Maxwell-Cattaneo type law is introduced in the system 
of hydrodynamic equations. Emphasis is made on the fact that the stabilization is only due to the 
difference in characteristic times for heat flux relaxation and instabilities onset. This gives further 
evidence that Eckart's like constitutive equations are responsible for the first order sytem exhibiting 
unphysical behavior. 
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I. INTRODUCTION 



The transport equations for the relativistic fluid have been well established within the 
framework of irreversible thermodynamics 

HQ 

. However, it was shown that Eckart's con- 
stitutive equation, as a closure equation for the heat flux including an acceleration term as 
driving force [3|, renders the system unstable |4| and seems to allow propagation of signals 
at speeds larger than the speed of light [5]. This instability and acausality in the relativistic 



gas lead to the formulation and use of higher order or extended theories [6H9|. 

The definite proof of the violation of Onsager's regression of fluctuations hypothesis was 
firstly shown by Hiscock and Lindblom eventhough they interpreted this contradiction with 
the tenets of irreversible thermodynamics as an early onset of instabilities [4] . In their paper 
it is shown that fluctuations in the linearized system of equations grow exponentially with 
a very small characteristic time. Meanwhile, buried deep in the fundamental hypothesis of 
non-equilibrium thermodynamics, Onsager's regression assumption assert that spontaneous 
fluctuations of microscopic origin should decay following the linearized equations for the 
state variables. It is this assumption that is violated when Eckart's equation for the heat 



flux is used, as it is shown in Ref. 



10]. 



As a result, Hiscock and Lindblom suggested in their work that, due to the unphysical 
behavior predicted by Eckart's formalism, "standard" first order theories should be discarded 
in favor of the second order one developed by Israel. Even though the causality in the 
extended theories was at the time verified, to the authors' knowledge the stability within 
such framework was not assessed. 

In this work we revisit the stability analisis for fluctuations within the linearized system 
by using a Maxwell-Cattaneo type constitutive equation for heat. This calculation not only 
validates the suggestion in Ref. J4| but also sheds light on the "stabilization" mechanism 
by showing that the dynamics of the fluctuations are dominated by two competing drives: 
Eckart's constitutive equation and Cattaneo's damping. 

The rest of this paper is divided as follows. In Sect. II the system of transport equations 
is established and the Maxwell-Cattaneo law for the heat flux is introduced. The equations 
are linearized for fluctuations around equilibrium values for the state variables in Sect. III. 
The dispersion relation and its cualitative analysis are shown in Sect. IV, while Sect. V 
includes the corresponding discussion and final remarks. 
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II. RELATIVISTIC TRANSPORT EQUATIONS 



The transport equations for a simple relativistic fluid are obtained from the conservation 
equations for the particle and momentum-energy fluxes, that is 

K» = o (i) 

and 

T£ = (2) 
where is the particle four flux given by 

iV M = (3) 

and 

TIF I I 

Tjf = —u»u v + pK + W u + -q»u v + -v?q v (4) 
c c c 

is the stress-energy tensor. Here n is the particle number density, u^ 1 the hydrodynamic four- 
velocity, c the speed of light, e the internal energy density per particle and p the hydrostatic 
pressure. The dissipative fluxes are the heat flux q v and the Navier tensor 11^. The spatial 
projector h£ for a + + +— signature is given by 

K = V + ^ ( 5 ) 

which satisfies h^u u = for // = 1,...,4, and the orthogonality conditions implied in this 
3+1 representation are 

UfM W u = 0, qX = (6) 
The introduction of Eq. ([3]) in Eq. ([1]) yields the relativistic continuity equation, i. e. 

n + n6 = (7) 

with 9 = u u u . Both energy and momentum balances are extracted from Eq. (j2J). For 
v = 1, 2, 3 one obtains the momentum balance 

(™ + 4) Uu + ( ™ + ~ 2 Q J U v +p,X + K;„ + \ {qfab, + q^ + 6q u + u^q^) = (8) 

Considering v — 4 leads to a total energy balance from which the mechanical energy needs 
to be substracted in order to establish the heat flux. A shortcut leading directly to the 
internal energy equation is to consider the projection u v T£ = which leads to 



ne + p6 + u^IIJJ + g£ + ^u v q v = (9) 

Using the relation e = (§^) T n + (J|;) T together with the ideal gas law p = nkT , which 
holds in the relativistic case [l|, one readily obtains an evolution equation for T, i. e. 

nC n f + p6 + u^ u + g£ + ^u v q v = (10) 

The set of hydrodynamic equations for the relativistic fluid is thus given by Eqs. ([?]), (jSJ) 
and (ITUj) . The set is not complete and constitutive equations have to be introduced in order 
to express the dissipative fluxes in terms of the state variables. In 1940, Eckart proposed, 
from a purely phenomenological approach by enforcing consistency with the second law 
of thermodynamics, a coupling of heat with hydrodynamic acceleration additional to the 
Fourier term [3] . The resulting system of equations was shown to yield unphysical behavior 



by Hiscock and Lindblom (4]. Following their suggestion, we here study the behavior of the 
linearized system by introducing for the heat flux, a Maxwell-Cattaneo constitutive relation. 
Such relation consists in an evolution equation for the heat flux, that is 

rq u + q u = -Kh^ {T^ + ^u)j (11) 

where r is a relaxation time such as r = Dt /CL with C s the speed of sound. Notice that 
Eq. (TTTj) becomes Eckart's constitutive equation [3] if r = 0. 



III. LINEARIZED TRANSPORT EQUATIONS 

In this section a set of linearized equations will be obtained by assuming small fluctuations 
of the state variables around their equilibrium values (denoted by a subscript 0), that is 

F = F + 5F (12) 

Once this substitution is made, second (and higher order) terms on fluctuations are neglected. 
Then, for the continuity equation we have 

Sh + n 56 = (13) 

and for the momentum balance 
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[n e + po) 5u v + kT5n v + nk5T u - (58 <v - 2-q (8t£), + 6q u = 



(14) 



where ( and rj are the bulk and shear viscosity respectively; while the heat equation in this 
linear approximation reads 



nC n 5T + n kT 58 + 5q v v = 
Finally, Maxwell-Cattaneo's equation in terms of fluctuations is given by 

T 



(15) 



T Sq u + 5q v = -nh^ ST V + -^Su 



(16) 



Notice that Eqs. (|13J) and (fT5|) only depend on the velocity through 58. Because of that, it is 
convenient to calculate the divergence of Eq. (Il4p . Following this standard procedure, the 
longitudinal velocity gradient fluctuation mode is uncoupled and the momentum balance is 
reduced to a scalar equation for 58 



-5 (n £o + Po) S9. v + kT V5n + n kV 2 5T - C + ^ V 56 







(17) 



c ' y 3 

Now both the internal energy and momentum linearized balance equations depend on the 
divergence of the heat flux fluctuations. Thus, if one calculates de divergence of the linearized 
Maxwell-Cattaneo equation, Eq. (116|) . the system is reduced to a system of 4 scalar equations 
for the unknowns 5n, 58, 5T and 5q u u . The cualitative analysis of such system is carried out 
in the next section in the Fourier-Laplace space in order to address the behavior of such 
fluctuations. 



IV. THE DISPERSION RELATION 

In order to establish the dispersion relation, the system matrix is formed with Laplace 
and Fourier transforms of equations (TIB"]) . ()15j) and (I17j) and the divergence of Eq. (jl6j) . that 
is: 



.4 



s tiq 

-kT q 2 n ms + q 2 (( + f ) -n kq 2 j, 
nokT C n n s 1 



sTqk 
r .2 



(18) 



-nq 2 ts + 1 J 
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with s and q the Laplace and Fourier parameters respectively. The corresponding dispersion 
relation for the system can thus be written as 



d^s 4 + d^s 3 + d 2 s 2 + dis + d = 
where the coefficients do, d\, d 2 y d^ are given by: 



(19) 



di 



kT 



ni 



1 + 



do 



kT n 4 

n 9 

C n mn 



q 2 + 



no 



4 



K 



3 / C n n 



d 2 = q 2 



C + ~v 

mn \ 3 



C n no 



2kT n 



+ 



kT n 



m 



k 



(20) 

(21) 

(22) 

(23) 

(24) 
e*mno 

In order to determine approximate solutions for the fourth order dispersion relation, we 
follow the same ideas as in Ref. jll|. Since d^ is a small quantity, three of the roots of Eq. 
(|T9|) are calculated by neglecting the fourth order term and solving for 



mn Q 



d\ 



T 



d 3 s 3 + d 2 s 2 + dis + do = 



(25) 



which, using Mountain's method 12|, can be shown to yield a real root given by 

D T 

Sl 7 

and the two imaginary solutions 

D v D T /5 2 \ 2 . [5kf 



(26) 



(27) 



where use has been made of the fact that for an ideal gas, for z « 1, C n = |/c, 7 = | and 
the definitions Dy = (C + -£) and Dt = — for the viscous and thermal diffusivities 
have been introduced. These roots give rise to the usual behavior of density fluctuations, 
that is, a decaying mode with a characteristic time given by s\ and two oscillating ones with 
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frequencies given by the imaginary part of S2 and S3. Notice that the only relativistic correc- 
tions to this behavior corresponds to the Stokes-Kirchoff coefficient, the term in parenthesis 
in Eq. (1271) . Clearly, since the characteristic time r is of the order of the mean collision time 
~ 10~ n s and typical thermal diffusivities are in the 10~ 5 — 10 m 2 s -1 range, the first term 
dominates and the relativistic correction yields only a slight increase in the real part of the 
conjugate roots which vanishes as z — » in the non-relativistic limit. It is also important 
to point out the fact that the relativistic heat diffusion term, zD?, has opposite sign and, 
although small, would actually yield a small growth in amplitude. 

To approximate the fourth root, we use the information of the previous solutions by 
assuming they are still approximate roots of the complete equation, which is plausible given 
the smallness of d^. Following such proposal, we use the fact that the sum of all roots in 
an n-th order polynomial is equal to the ratio of the coefficients of the n — 1 and the n-th 
power to obtain 



3D T z\ 1 



S4 = - [T 2 c 2 J 



1 + -^—4 DyT^K + - -2— - (28) 



: + -c 2 T kt + — + -c A kn T T 2 

O O O III o 

Notice that, by taking r = 0, the fourth root is positive and thus corresponds to a growing 
mode. In Ref. [IX | . the existence of this positive large solution was used to argue against 
Eckart's constitutive equation, in particular the coupling of heat with acceleration. For 
t 7^ 0, we can approximate 

S4 ~-- + 9 2 (i>r+~^M (29) 
t \ 6 m J 

Now, Eq. f[2"9|) sheds more light on the issue since, by using the fact that r ^> one 
can conclude not only that the system of hydrodynamic equations with Cattaneo's closure 
(Eq. f)16p ) is stable but also realize that this stabilization takes place as the relaxation term 
dominates over the exponential growth caused by Eckart's coupling. In other words, the time 
derivative term in Maxwell-Cattaneo's constitutive equations yields not only a modified time 
scale for the evolution of the fluctuations but, more importantly, changes the overall physical 



behavior because o 



root found in Ref. ll| 



the different signs of the terms in the denominator of Eq. ( 129]) . The 



c 4 mn 

S4 7fr~ I 30 ) 
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yields an exponential growth in the structure factor while the result in the current calculation 

s 4 ~ -- (31) 

T 

corresponds to a finite spectrum which could eventually be observed and measured. 



V. SUMMARY AND FINAL REMARKS 

In the previous section, a dispersion relation for the system of linearized hydrodynamic 
equations for the relativistic fluid was obtained with a Maxwell-Cattaneo heat evolution 
equation. Such procedure produced a fourth order equation for the Fourier variable s for 
a given wavenumber q. This result is similar to the one obtained in Ref. but only in 
structure since, as was emphatically pointed out in Sect. IV, the overall behavior of the 
system changes radically. This can be clearly seen by comparing the approximate roots in 
Eqs. (I30p and (}3"Tj) . The former yields an exponential growth in the corresponding structure 
factor which destroys, theoretically, the spectrum which is clearly unacceptable. However, 
the latter predicts a decay in time of the fluctuations with a characteristic time determined 
by r, the relaxation parameter introduced by Cattaneo. 

In this way, the stabilization of the system predicted intuitively by Hiscock and Lindblom 
in 1985 is confirmed. To our knowledge, only the causality issue had been addressed so far. 
And moreover, the mechanism for the stabilization can be qualitatively attributed to the 
competing terms in the denominator of s 4 as given by Eq. (129"j) . By confirming that the 
use of a relaxation equation for the heat flux eliminates the instability in the system of 
equation by dominating over the runaway solution that is introduced by the acceleration- 
heat coupling, verifies the fact that it is such relation the one which in the first place led 
to the unphysical behavior of thermodynamic fluctuations in the relativistic gas. We thus 
conclude that the stability of the fluid within Maxwell-Cattaneo's formalism is in complete 
agreement with the claim that heat cannot be coupled with acceleration. 
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